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ABSTRACT 

The detection of mixed modes in subgiants and red giants by the CoRoT and Kepler space-bome missions allows us to investigate 
the internal structure of evolved low-mass stars, from the end of the main sequence to the central helium-burning phase. In particular, 
the measurement of the mean core rotation rate as a function of the evolution places stringent constraints on the physical mechanisms 
responsible for the angular momentum redistribution in stars. It showed that the current stellar evolution codes including the modelling 
of rotation fail to reproduce the observations. An additional physical process that efficiently extracts angular momentum from the core 
is thus necessary. 

Our aim is to assess the ability of mixed modes to do this. To this end, we developed a formalism that provides a modelling of the wave 
fluxes in both the mean angular momentum and the mean energy equations in a companion paper. In this article, mode amplitudes are 
modelled based on recent asteroseismic observations, and a quantitative estimate of the angular momentum transfer is obtained. This 
is performed for a benchmark model of 1.3 Mq at three evolutionary stages, representative of the evolved pulsating stars observed by 
CoRoT and Kepler. 

We show that mixed modes extract angular momentum from the innermost regions of subgiants and red giants. However, this transport 
of angular momentum from the core is unlikely to counterbalance the effect of the core contraction in subgiants and early red giants. In 
contrast, for more evolved red giants, mixed modes are found efficient enough to balance and exceed the effect of the core contraction, 
in particular in the hydrogen-burning shell. Our results thus indicate that mixed modes are a promising candidate to explain the 
observed spin-down of the core of evolved red giants, but that an other mechanism is to be invoked for subgiants and early red giants. 
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1. Introduction 


The Co RoT (|Baglin et al. |2006a|b[ [Michel et al.||20()8) l and 


Kepler (Borucki et al. 20l6|l space-borne missions provided a 


wealth of observed stars exhibiting solar-like oscillations, from 
the main-sequence to the red-giant phases (see jChaplin & Miglioj 
[2013 for a review). A large number of those stars are low-mass 


evolved stars showing a rich spectrum of mixed modes, which 
behave as acoustic modes in the stellar envelope and as grav¬ 
ity modes in the core. They are therefore detectable at the stel- 


lar surface while yielding information on the innermost regions 

(e.g., Dziembowski|1971 

Scuflaire|1974[ Aizenman et al.|1977 

Dziembowski et al.(2001 

Dupret et al.|2009[). 


While mixed modes allowed unveiling the structure of the 
red giants core ([Bedding et al.|20lT] Mosser et al.|2011 [2014 1 , 


the detection of associated rotational splittings enabled mea- 
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suring the mean core rotation of subgiant and red giant stars 


([Beck et al.[[2012 Deheuvels et al.p012[ [Mosser et al.[[2012bl 


[Deheuvels et al. 2014| l. It turns out that the core of low 
red giant stars slows down during evolution. Since local con¬ 
servation of angular momentum would imply a spin-up of the 
core, these observations indicate a strong transfer of angular 
momentum from the inner to the outer layers of these stars. 
However, current stellar evolutionary rotating models that take 
the transport of angular momentum into account fail to repro¬ 
duce the observations. Meridional circulation or shear instabil¬ 
ities are not efficient enough to slow down the red-giant cores 
(e.g., [Eggenberger et al.| 20T^ Marques et al. 2013 Ceillier 
[et al.[(^13| l, and recent attempts to include the effect of prop¬ 
agative internal gravity waves ( jEuller et al.|2014)l or magnetic 
fields using the Taylor-Spruit dynamo formalism (Cantiello et al. 
[2014[ ) failed to solve the problem. Nonetheless, for subgiants 
and early red giants, Rudiger et al. (2015 1 have recently shown 


that magneto-rotational instabilities of a toroidal magnetic field 
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could explain the angular momentum redistribution in subgiants 
and early red giants (see also Maeder & Meynet|2014| l. 

Hence, an additional physical mechanism that efficiently ex¬ 
tracts angular momentum from the core of red giants seems to be 
needed. Non-radial modes in subgiants and red giants are poten¬ 
tial candidates. Indeed, prograde and retrograde mixed modes 
are differentially damped in the presence of rotation, allowing 
a net transport of angular momentum (e.g., |Ando|[T98^ |Lee &| 
Saio||1993 1. Our aim is thus to investigate the ability of mixed 


modes to transfer angular momentum from the inner radiative 
interior to the outer convective layers. 

In the hrst article of this series (hereafter Paper I), we used 
the transformed Eulerian mean (TEM) formalism to account for 
the effect of wave driving through both the wave momentum flux 
in the mean angular momentum equation and the wave heat flux 
in the mean energy equation. The wave held was then modelled 
using the asymptotic, quasi-adiabatic, and slow-rotation approx¬ 
imations, valid in the dense radiative layers of evolved stars. This 
allowed us to obtain an explicit expression of the wave huxes ap¬ 
pearing in the mean angular momentum and energy equations. 
This article, the second of this series, is dedicated to a quanti¬ 
tative determination of the angular momentum transported by 
mixed modes, which requires the knowledge of mode ampli¬ 
tudes. Recent CoRoT and Kepler observations ( |Baudin et al.| 
|201 1[ [Mosser et al.|2012^ [Samadi et al.||2012) l allow us to de¬ 
termine mixed-mode amplitudes and subsequently to determine 
the efficiency of angular momentum transport by mixed modes. 

The paper is organized as follows: Sect, [^recalls the formal¬ 
ism developed in Paper I, while Sect, [^presents the modelling of 
the equilibrium models and related oscillations. In Sect.|^ based 
on the observations, we describe the modelling of mode ampli¬ 
tudes. Finally, the estimate of the angular momentum transport 
efficiency is presented in Sect.|^ and Sect.j^is dedicated to dis¬ 
cussions and conclusions. 


2. Theoretical formalism 

In this section, we recall the formalism developed in Paper I as 
well as the related assumptions for both the mean flow and the 
wave fluxes appearing in the mean equations. 


where ur = r sin 0, the overbar is the azimuthal average, is the 
residual meridional velocity, Vy is the radial eddy viscosity, s is 
the specihc entropy, and the prime denotes perturbations associ¬ 
ated with the non-radial oscillations, so that v^, v', v'g are the az¬ 
imuthal, radial, and meridional component of the wave velocity 
held, and s' the wave Eulerian perturbation of entropy. We have 
also introduced the Lagrangian derivative d/df = d/dt + r djdr. 

Under the same assumptions as for Eq. ([T]), the mean entropy 
equation is 


<P) 


d_W 

df 


yl Qy 


{r^S) + {Q) 


( 5 ) 


with 


and 


= <p) i' v;, 




(6) 


( 7 ) 


where T is the temperature, e is the nuclear energy generation 
rate and;^f is the thermal conductivity (see, e.g.. Sect. 6 of Mathis 
& Zahn|f2004 and also Maeder & Zahn|[i998 1. This expresses 
energy conservation on level surfaces. 

To quantify the effect of the waves on the rotation prohle, we 
further need to provide a modelling of the wave heldj 


2.2. Wave flux modelling 

Our aim is to determine the angular momentum redistribution 
induced by mixed modes in the radiative region of low-mass 
evolved stars. Therefore, the wave Huxes in Eqs. Q and Q must 
be modelled. This has been done in Paper I, so that for the mo¬ 
mentum equation we have 
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2.1. Mean flow equations 

Following Paper I, we adopt the transformed Eulerian mean for¬ 
malism (TEM). In addition, we assume shellular rotation (see 
|Maeder|[2009| for an extensive discussion) so that at dominant 
order the mean angular momentum equation is 


d(r2Q„)_ 

“ y2 Qy ['■ 


dt 


: + ^Tsh ear + ^waves)], (1) 


where Qq is the rotation angular frequency that depends only 
on the radius, p is the density, the overbar and the symbol () 
denotes the horizontal average {i.e. azimuthal and meridional). 
The Huxes are deHned by 


Tcirc = - ^ <p>r-Qo Uj , 

^hear “ \P) rv » 

or 


Tw; 


<P){CT 


v^ Vj -H 2 cos 0Qo VgS' I 


( 2 ) 

( 3 ) 

( 4 ) 


where the coefficients S™, C”', are given by Eqs. (A.25) 
to (A.28) of Paper I, and o’ = cr^ H- mQo (with ctr the modal 
frequency). The amplitude corresponds to the amplitude of 
a mode of a given angular degree, £, and azimuthal order, m. 
We note that the amplitude is considered statistically constant in 
time since we deal with solar-like oscillations that result from a 
balance between mode driving and damping (see Samadi|201 1 
for details). 

For the entropy equation, we obtained 




( 9 ) 


where 


Anr^pT \ V 

,2 ^(^+1) 


/V^_,WcU 




(10) 


' In this article, by wave we denote global mixed modes. 
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Fig. 1. Evolutionary track on the Hertzsprung-Russell (HR) di¬ 
agram of a 1 3Mq model, showing the location of the selected 
models (see Sect.j^for details). 


with L the luminosity, N the buoyancy frequency, V and Vad the 
actual and adiabatic temperature gradients, respectively, and 
the radial component of the eigenfunction for a mode of angu¬ 
lar degree £ and azimuthal order m. The overbars and angular 
brackets have been dropped for ease of notation. 

We also recall that in the course of deriving Eqs. ([8])and0 
several simplifying approximations were performed, namely: 


Table 1. Characteristics of the selected models as described in 
Sect. 1^ Model MO is a subgiant, model Ml is at the base of 
the red-giant branch, and model M2 is a red giant. Here, Av is 
the large separation and Vmax is the frequency at the maximum 
height of the oscillation power spectrum. 


Model 

Age (Myr) 

R/Re 

L/Lq 

Av (pHz) 

l^max O^Hz) 

MO 

4401 

2.47 

4.95 

39.5 

666 

Ml 

4590 

3.55 

6.78 

22.9 

339 

M2 

4791 

9.51 

38.3 

5.24 

48.8 



r/R 


Fig. 2. Synthetic rotation profiles, computed as described in 
Sect. |3.2| for models MO, Ml, and M2 as a function of the radius 
normalised by the total radius. We used logarithmic scales. 


1. The quasi-adiabatic approach: It consists of neglecting the 
difference between adiabatic and non-adiabatic eigenfunc¬ 
tions in the full wave equations. This approximation is valid 
provided the local thermal timescale is much longer than 
the modal period. This is justified in the radiative region of 
evolved low-mass stars. 

2. The low-rotation limit: We consider the low-rotation limit in 
the radiative regions or, more precisely, ctr » Qq- This is 
justified by recent inferences of the rotation rate in the core 
of subgiants (|Deheuvels et al. 2012 20141 and red giants 
(Mosser et al. 2012b| l using seismic constraints from Kepler. 

3. The asymptotic limit: We used an asymptotic description for 
gravity modes (e.g., |Dziembowski et al.|2001[ Godart et al. 
|2009^ , which is valid for mixed modes in the inner radiative 
region of subgiants and red giants (e.g., Goupil et al.|2013|l. 


Consequently, our formalism is valid in the dense radiative core 
of evolved slowly rotating stars. 


3.1. Equilibrium models and related eigenfunctions and 
eigenfrequencies 


We used the stellar evolution code CESTAM (Marques et al. 
20131 to compute the equilibrium models. The atmosphere 
was computed assuming a grey Eddington approximation. 
Convection was included according to |Canuto et al.| ( [T996| l, with 
a mixing-length parameter a - 0.67. The initial chemical com¬ 
position follows Asplund et al. (2005 1 , with a helium mass frac¬ 
tion of T = 0.261 and a metallicity of Z = 0.0138. We used 
the OPAL equation of sta te ( Rogers et al.||1996 l and opacities 
(Iglesias & Rogers 1996 1 , complemented, at T < 10"^ K, by 
the [Alexander & Fergusoti] ( |1994| l opacities. We adopted the 
NACRE nuclear reaction rates from Angulo et al.| ( 1999 1 except 
for the -H p reaction, where we used the reaction rates given 
in |Imbriani et al.| ( |2004| l. 

Finally, we used the ADIPLS code to compute adiabatic os¬ 
cillations ( Christensen-Dalsgaard|20()8l |201 l| l, that is, to obtain 
the eigenfunctions and eigenfrequencies. 


3. Computing selected stellar models 

We have computed an evolutionary sequence of a benchmark 
model of a 1.3 Mq evolved low-mass star representative of 
CoRoT and Kepler observations. We selected three M - 1.3Mo 
models at three post-main sequence stages of this evolutionary 
sequence, as discussed below. The location of the models in the 
Hertzsprung-Russell (HR) diagram is shown in Fig.[T] and their 
characteristics are summarised in Table [T] 


3.2. Rotation profiles 

We consider that rotation has negligible effects on the equilib¬ 
rium structure of the star and that it can be considered as a per¬ 
turbation for the oscillations. This approximation is consistent 
with the formalism developed in Paper I. 

The rotation profiles are therefore determined by computing 
synthetic profiles that match the observations. The synthetic pro¬ 
files should have a core that rotates faster than the envelope, with 
a transition region located at the H-burning shell in the case of 
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Table 2. Parameters of the synthetic profiles as described in 


Sect. |3.2| (see also Eq. 111. 


Model 

Plclln (nHz) 

Qs /27r (nHz) 

rtIKav 

w//?star 

MO 

600 

200 

0.1452 

1.452 10-2 

Ml 

1000 

55 

1.009 10-2 

1.211 10-2 

M2 

500 

8.1 

3.316 10-2 

3.015 10-“ 


red giant stars (models Ml and M2). The width of the transition 
region is related to the width of the H-burning shell. We conse¬ 
quently used the following profile 


Q(r) — Q5 -H 



1 -H erf 



( 11 ) 


where Qs and Qc the surface and central angular velocities, 
respectively, and r, and w are the radius and width of the burning 
shell. Qc was chosen to reproduce the mean core rotation rate of 
evolved low-mass stars as derived from Kepler observations (see 
[Mosser et al.|2012b{ [Deheuvels et al.|2014| l, and the surface ro¬ 
tation rate (Qs) was obtained so as to conserve the total angular 
momentum from the subgiant stage (model MO). 

For the subgiant model (model MO), we also used Eq. 
but with both r, and w given by the radius of the central convec¬ 
tive zone at the end of the main sequence (just before it disap¬ 
pears). In this case, the central and surface rotation (Qc and Qs, 
respectively) were chosen to reproduce the results of|Deheuvels 


et al. (|2012l. The values of the parameters for models MO, Ml, 


and M2 are summarised in Table and the resulting rotation 
profiles are displayed in Fig. 


4. Modelling mode amplitudes 



V [/^Hz] 



200 250 300 350 400 450 500 

V [/.iHz] 


There are mainly two approaches to compute mode amplitudes. 
The first is based on a full non-adiabatic computation includ¬ 
ing a time-dependent treatment of convection. This procedure is 
time-consuming, however, and still suffers from theoretical un¬ 
certainties (see Dupret et al.|2009t Grosjean et al.|20f4 for de¬ 
tails). The second is based on recent CoRoT and Kepler obser¬ 
vations that allowed us to measure the amplitudes of solar-like 
modes in a wide variety of stars in different evolutionary stages. 
This wealth of observations made it possible to establish scaling 
relations that provide mode amplitudes versus global stellar pa¬ 
rameters (e.g., [Mosser et al.|2012a[ Samadi et al.|2012l l. In this 
paper, we adopt the latter approach. 


4. 1. Radial mode amplitudes 

We first express the amplitudes of radial modes. The surface ve¬ 
locity of the radial modes, at the frequency of maximum height 
(Vmax) in the oscillation power spectrum, follows the scaling re¬ 
lation (see jBelkacem & Samadi|2013| for details) 

Vr=0 (V = Vmax) _ / Teff / Vmax / Av 

ye=o,o (y - ^max) \^G,eff/ \ 1^0,max / \Avq/ 

( 12 ) 



20 30 40 50 60 70 80 

u [/tHz] 


Fig. 3. Mode amplitudes versus mode frequencies for models 
MO (top panel). Ml (middle panel), and M2 (bottom panel). The 
amplitudes are computed as described in Sect. 


where Tes is the effective temperature, Av is the large separation, 
and the O symbols stand for the solar values. Here we adopt the 
values To^eff = 5777 K, Vo_niax = 3050jt/Hz, and Avq = 135yuHz. 

To account for the frequency dependence of radial mode am¬ 
plitudes, we follow the work ofjMosser et al.|(|2010|l , who used 


a Gaussian envelope centred on Vmax with a full-width at half¬ 
maximum dVenv. that is, 

via M = via (y = ^max) ^-(-’'■””>4(2-4 , (13) 
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with cr = 5venv/(2 V2 111 2). Based on Kepler observations, we 
take dVenv = 150/tHz for MO, dVenv = 70/rHz for Ml, and 
5venv = 30 /rHz for M2 for the following. 


4.2. Non-radial mode amplitudes 

The next step is to consider the amplitudes of non-radial modes. 
The ratio between radial and non-radial mode amplitudes of 
comparable frequencies can be written as (see |Benomar et al. 
|20T4l for details) 








(=0 

r7 


(14) 


where A4f - f ^^dm/^^(f?) is the mode mass, ^ is the eigen¬ 
function, m is the mass profile, R is the total radius, and Tf is the 
mode linewidth that is related to the mode damping rate rjt by 
Vi = rjdn. 

For high angular degrees, the radiative work of the oscilla¬ 
tions in the innermost regions dominantes the work done in the 
uppermost layers (see Grosjean et al.|[20T4| for details). It fol¬ 
lows that the mode amplitudes of high angular degrees become 
very small due to strong radiative damping. We thus dehne a cut¬ 
off angular degree, fmax, below which we assume that the work 
of non-radial modes equals the work of radial modes 


r^=oA1f=o - 


t<tm 


(15) 


where the product between the mode mass and the mode damp 
ing rate equals the work integral (e.g., Dupret et al. [2009 


Benomar et al.||2014[[Grosjean et al.|2014|l. Now, using Eq. 

together with Eq. (|14f leadsTo 


© 


1/2 

‘-max 


(v) Me= 


(y) M{<t^ 


(16) 


For { > ^max, we assume that the radiative damping in the inner¬ 
most region becomes dominant over the damping of the upper 
layers so that mode amplitudes become negligible. 

The cut-off angular degree is obtained by assuming that the 
work done over an oscillation period through the effect of ra¬ 
diative losses in the innermost layers becomes of the same or¬ 
der of magnitude as the work done in the uppermost layers. We 
obtain - 6 for MO, fmax = 4 for Ml, and - 2 for 
M2 (see Appendix for details). Note that this approach for 
the modelling of mode amplitudes provides a lower limit to the 
amount of angular momentum transported by mixed modes since 
we neglect the contribution of angular degrees higher than 
However, this is enough to determine whether mixed modes are 
able to slow down the rotating core of red giants. 

Finally, we assume that the relation between the surface 
mode velocity and the mode amplitude is given by = 
a? cr^l^^P/2 (see Samadi 20111. The results obtained for mode 


amplitudes of angular degrees t - {0,1,2) are presented in Fig.[^ 
for the three models MO, Ml, and M2, described in Sect.j^ 


5. Quantitative estimate of anguiar momentum 
extracted by mixed modes 

Our aim is to estimate the rate of angular momentum transported 
by mixed modes for the benchmark models MO, Ml, and M2. To 
this end, as in Paper I, we introduce the notation 

1 

7 = -^-(r27-waves) , (17) 



r/Rc 


Fig. 4. Rate of temporal variation of the mean specific angular 
momentum (7 as dehned by Eq. 17 1 as a function of the star 
radius normalised by the radius of the base of the convective 
envelope. Red corresponds to 7 < 0 and blue to 7 > 0. 


where the right-hand side is given by Eq. (|^. Using the approach 
presented in Sect. to derive the amplitudes, we compute the 
wave fluxes in the right-hand side of Eqs. Q and 0 using 
Eqs. (|^ and Q. 

The contribution of the mixed modes to the mean angular 
momentum equation is shown in Fig.|^for models MO, Ml, and 
M2. The main effect of mixed modes is to slow down the core 
of the stars in the region near the maximum of the buoyancy fre¬ 
quency that also corresponds to the edge of the fast rotating core. 
For the red-giant models, this region is the hydrogen-burning 
shell. The decrease shown towards the upper radiative layers is 
the result of both a lower rotation rate and a lower value of the 
buoyancy frequency. This behaviour is similar in the three mod¬ 
els. In addition, the amount of angular momentum extracted by 
mixed modes hardly depends on the rotation gradient between 
the core and the envelope (the sensitivity of the results to the 
parameter w in Eq. GD remains weak). 

The amount of angular momentum extracted by mixed 
modes increases with the evolutionary stage of the star. This ef¬ 
fect is the result of several factors. First, as shown in Fig.[^ mode 
amplitudes increase from models MO to M2 and thus more en¬ 
ergy is available to transport angular momentum. Moreover, the 
buoyancy frequency also increases and so does the radial wave 
number in the gravity-mode cavity. As the energy exchanges be¬ 
tween modes and the background is proportional to the radial 
wave number, the amount of angular momentum extracted also 
increases. Finally, the number of mixed modes between two p- 
dominated modes significantly increases between MO and M2. 

Furthermore, the angular momentum is mainly carried away 
by g-dominated modes of the highest angular degrees. Indeed, 
these modes are very efficiently trapped in the core and undergo 
a larger radiative damping (due to the larger radial wave num¬ 
ber). The predominance of high angular degree modes implies 
that 7 is sensitive to our determination of 7max- As mentioned 
before, the estimate of the amount of angular momentum ex¬ 
tracted in the radiative regions is thus a lower limit. In other 
words, adding the contribution of angular degrees higher than 
^max would Strengthen these trends. 





























6 


Belkacem: Angular momentum redistribution by mixed modes in evolved low-mass stars 


To proceed and assess the efficiency of the transport of an¬ 
gular momentum for our benchmark models, we compare two 
timescales, namely: 


- the timescale associated with the efficiency of the transport 
of angular momentum by mixed modes, defined as 


j 

pr^D.0 


(18) 


- the timescale associated with the contraction of the star, de¬ 
fined as 


1 d 
dr 


{pr‘^Q.0 f) 



(19) 


where r - dr/df. 

The two timescales are shown in Fig.|^(top panel). In all the 
radiative regions of models MO and Ml, the extraction of angular 
momentum remains negligible compared to the star contraction. 
In contrast, for model M2, the timescale of angular momentum 
extraction is of the same order of magnitude as the contraction 
timescale and even lower in the hydrogen-shell-burning region, 
for which the contraction is maximal (see Fig.|^ bottom panel). 
However, in the very centre or in the upper radiative region, the 
mixed mode extraction of angular momentum is again negligi¬ 
ble. 

Therefore, we conclude that the extraction of angular mo¬ 
mentum by mixed modes is negligible in subgiants and early red 
giants, but becomes strong in the hydrogen-burning shell in stars 
higher on the red-giant branch. In these cases, mixed modes are 
able to counterbalance the spin-up due to the star contraction 
and can thus enforce a spin-down in those layers. We note, how¬ 
ever, that the exact location in the HR diagram of the transition 
between inefficient and efficient extraction by mixed modes is 
likely to depend on the mass and internal physics of the models. 


6. Concluding remarks 


Based on the theoretical formalism developed in paper I, we 
computed the amount of angular momentum extracted by mixed 
modes in evolved low-mass stars. To this end, we considered 
1 3Mq benchmark models of a subgiant (model MO), an early red 
giant (model Ml), and a more evolved red giant (model M2) for 
which we computed realistic mode amplitudes based on recent 
CoRoT and Kepler observations. We found that mixed modes 
extract angular momentum in the innermost region of subgiants 
and red giants. This extraction of angular momentum is found 
negligible for the subgiant and lower RGB models, whereas it is 
efficient enough to counterbalance the effect of the stellar con¬ 
traction for the most evolved model. 

For this evolved star, mixed modes efficiently transport the 
angular momentum in the hydrogen-burning shell (that corre¬ 
sponds to the maximum of the buoyancy frequency) but becomes 
less efficient in the very central layers. We note, however, that 
mixed modes can be considered as an external torque for merid¬ 
ional circulation (e.g., Mathis et al.|2013| . Therefore, such a lo¬ 
calized spin-down in the hydrogen-burning shell is likely to feed 
a strong redistribution of angular momentum in the central lay¬ 
ers by meridional circulation. This would nevertheless deserve a 
quantitative estimate. 

Indeed, evolutionary calculations including the effect of 
mixed modes coupled with meridional circulation and shear in¬ 
stabilities are the next step for a definitive conclusion on the abil¬ 
ity of mixed modes to slow down the cores of evolved red giants. 



r/Rc 



r/Rc 


Fig. 5. Top panel: Timescales versus normalised radius {i.e., 
normalised by the radius of the base of the convective envelope) 
for models MO, Ml, and M2. The solid lines correspond to the 
timescale associated with the transport of angular momentum by 
mixed modes (see Eq. 18 1 and the dotted lines correspond to the 
timescale associated with the contraction of the star (see Eq.[T^. 
Bottom panel: Same as the top panel, except that only a zoom 
for model M2 is shown. 


Modelling of mode amplitudes would require further investiga¬ 
tion since we assumed that the contribution of high-angular de¬ 
gree modes would be negligible due to their high radiative damp¬ 
ing. Therefore, we believe that our calculations provide a lower 
limit. Full non-adiabatic computations would provide more ac¬ 
curate results, however. 


We conclude that, for subgiants and early red giants, our re¬ 
sults indicate that extraction of angular momentum by mixed 
modes cannot explain the current observations ( [Deheuvels et al. 
20141. Other physical mechanisms such as internal gravity 


( |Talon & Charbonnel||2008[ [Fuller et al.||2014| l or mag¬ 
netic fields ( [Rudiger et al.|2015| ) ffiould also be present in these 
stars in order to extract enough angular momentum from the 
core. In contrast, as demonstrated in this work, the transport of 
angular momentum by mixed modes is a serious explanation for 
the observed spin down of the core of evolved stars on the red- 
giant branch. At least, given the order of magnitude of the in- 
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duced angular momentum extraction, it seems difficult to neglect 
this mechanism in the future. 

Finally, we note that an extension of this work to central he¬ 
lium burning stars would be desirable in the future. Indeed, even 
if the mode amplitudes are similar for stars on the vertical branch 
and on the clump for the same v^ax (e.g., Mosser et al.|2012a) l, 
the maximum of buoyancy frequency is lower and shifted to¬ 
ward higher radius for clump stars (e.g., Montalban et al.|2013| l. 
Therefore, we can expect a lower efficiency of the transport by 
mixed modes but this needs to be definitively quantified. 
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Appendix A: Cut-off angular degrees 


The cut-off angular degree, ^max as introduced in Sect. |4.2[ can 
be defined as the angular degree for which the work performed 
over an oscillation period through the effect of radiative losses in 
the inner-most layers becomes of the same order of magnitude as 
for the work performed in the uppermost layers. More precisely, 
^max is obtained at v = Vmax through the relation 

(^^inax^^4iax)upper layers (^^max^^^max)inner layers ‘ (-^-1) 

The left-hand side of Eq. (jATj can be identified with a radial 
mode of comparable frequency for which the mode mass is com- 



0, V = Vmax) = 0.95/rHz for Ml, and rjit - 0, v = Vmax) = 
0.06 for M2. 

The right-hand side of Eq. ( |A.l| l can be identified with a 
gravity dominated mode for which the mode mass is computed 
following Sect. and the damping rate is obtained using the 
asymptotic relation (e.g. [Dziembowski et al.|2001||Godatt et al. 
|2^ 


^ (^max? v) — 


[W(Cax + 

//‘k,dr 




f 


Vad - V VadA^gT 

V Pr^ 


dr. 


(A.2) 

where F is the pressure, L is the luminosity, is the base of the 
convective region, ctr - Inv, V is the temperature gradient and 
Vad its adiabatic counter-part. 

Einally, using Eqs. (A.l i, and (A.2 1 , one obtains ^max = 6 for 
MO, ^max = 4 for Ml, and fmax = 2 for M2. 



































